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The quarkonium states in a quark-gluon plasma is examined with the heavy quark potential 
implied by the holographic principle. Both the vanila AdS-Schwarzschild metric and the one with 
an infrared cutoff are considered. The dissociation temperature is calculated by solving the Schro 
dinger equation of the potential model. In the case of the AdS-Schwarzschild metric with a IR 
cutoff, the dissociation temperatures for J/tp and T with the U-ansatz of the potential are found to 
agree with the lattice results within a factor of two. 
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I. INTRODUCTION 

AdS/CFT duality opens a new avenue towards a qualitative or even semi-quantitative understanding of the nonper- 
turbative aspect of a quantum field theory [ij Q . Because of the isomorphism between the conformal group in four 
dimensions and the isometry group of AdS^ space, it is conjectured that the a string theory formulated in AdS^ x 
is dual to the = 4 supersymmetric Yang-Mills theory (SUSY YM) on the boundary. The latter is believed to 
be conformal at quantum level and the global 5t/(4) symmetry of its i?-charges corresponds to the isometry of S5. 
In particular, the low energy limit of the classical string theory, the supergravity in AdS^ x corresponds to the 
supersymmetric Yang-Mills theory at large Nc and large 't Hooft coupling. 

A S iVcff^M (1) 

Various field theoretic correlation functions can be extracted from the metric fluctuations of the gravity dual and the 
expectation value of a Wilson loop operator is related to the minimum area the loop spans in the AdS^ bulk. 

Notable success has been made in the application of the AdS/CFT duality to the physics of quark-gluon plasma 
(QGP) created by RHIC, even though the underlying dynamics of QCD is very different from that of a supersymmetric 
Yang-Mills theory [1,13, [1]. Among them are the viscosity-entropy ratio @, the jet quenching parameter [1| which 
are closer to the observed values than the perturbative results. While it is premature to conclude that every aspect 
of RHIC physics can be explained in terms of SUSY YM, the AdS/CFT duality provides unprecedented references 
since this is the only case where the strong coupling properties of a quantum field theory can be calculated reliably. 

The heavy quarkonium dissociation is an important signal of the formation of QGP in RHIC. In the deconfinement 
phase of QCD, the range of the binding potential between a quark and an anti-quark is limited by the screening length 
in a hot and dense medium, which decreases with an increasing temperature. Beyond the dissociate temperature, Td, 
the range of the potential is too short to hold a bound state and the heavy meson will melt. The lattice simulation of the 
quark-antiquark potential and the spectral density of hadronic correlators yield consistent picture of the quarkonium 
dissociation and the numerical values Td. It is the object of this paper to calculate Td using the heavy quark potential 
of iV = 4 SUSY YM extracted from its gravity dual @ 0] 0- Although the potential model applies only in the 
non-relativistic limit which is not the case when the 't Hooft coupling of the = 4 SUSY YM, A becomes too strong, 
it can be justified within the lower side of the range of A used in the literature to compare AdS/CFT with the RHIC 
phenomenology, i.e. 

5.5 < A < 6tt. (2) 
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The upper edge is obtained by substituting into Q Nc = 3 and the QCD value of gyu at RHIC energy scale 
(^ym/I^'"') — 1/2) and the lower edge is based on a comparison between the heavy quark potential from lattice 
simulation with that from AdS / CFT ^ . 

We model the quarkonium, J/ijj and T as a non-relativistic bound state of a quark and an antiquark. The wave 
function for their relative motion satisfies the Schro dinger equation 

- TT^VV + V,s.{r)i^ = -Ei: (3) 
Im 

where m = A//2 is the reduced mass with M the mass of the heavy quark and E{> 0) is the binding energy of the 
bound state. Because of the screening of QGP, the effective potential energy has a finite range and is temperature 
dependent. The dissociation temperature of a particular state is the temperature at which its E is elevated to zero. 

In the next section we shall calculate the dissociation temperature semi-analytically with the heavy quark potential 
extracted from the vanila AdS-Schwarzschild metric. The case with an infrared cutoff will be examined in the section 
III and the section IV will conclude the paper. 



II. THE HEAVY QUARK POTENTIAL FROM ADS/CFT 



The free energy of a static pair of qq separated by a distance r at temperature T is given 

-iF(.,T) ^ tr < W^{L+)W{L^) > 

tr< W''{L+) >< W{L-) > ^ ' 

where L± stands for the Wilson line running in Euclidean time direction at spatial coordinates (0,0, ±ir) and is 
closed by the periodicity. We have 

W{L±) EE Pe-'' dtAo(tfi..O,±^r) (5) 

with Aq the temporal component of the gauge potential subject to the periodic boundary condition, Ao{t ^^r) — 
Ao(t, r). The trace in (O is over the color indexes and < ... > denotes the thermal average. The symbol P enforces 
the path ordering. The corresponding internal energy reads 

U{r^T)==-T^'4-(^^\ (6) 



dT \ T 

In the de-confined phase, the numerator of Q factorizes at large separation, i.e. 

lim < W\L+)W{L^) >=< W\L+) >< W{L^) > (7) 

r—^oo 

and therefore limr ^ oo F{r,T) = 0. The thermal average of a single Wilson line, < W{L) >, is independent of the 
spatial coordinates. 

Two ansatz of the potential model have been explored in the literature [lO|: the P-ansatz which identifies Vcs. with 
F{r,T) and the ?7-ansatz which identifies VcS. with U{r,T). The lattice QCD simulation reveals that the U ansatz 
produces a deeper potential well and thereby higher because the entropy contribution is subtracted. This remains 
the case with holographic potential as we shall see. 

According to the holographic principle, the thermal average of a Wilson loop operator 

iy(C) = Pe^'fc'^^''^"^^") (8) 

in 4D = 4 SUSY YM at large Nc and large 't Hooft coupling corresponds to the minimum area ^min. [C] of the 
string world sheet in the 5D AdS-Schwarzschild metric with a Euclidean signature, 

ds^ = n^T^y^fdt^ + dx^) + ^dy\ (9) 

r/ 

bounded by the loop C at the boundary, y — s- cx), where f — \ — . Specifically, we have 

tr < W{C) >= e-^^"""['^l. (10) 



3 



For the numerator of ([4]), C consists of two parallel temporal lines (t,0, 0,±^) and the string world sheet can be 
parameterized by t and y with the ansatz xi = X2 — ^ and 0:3 a function of y. The induced world sheet metric reads 

ds^ . .^T-y-fdt- + [-'T-v\^)\ ^^W- (11) 
Minimizing the world sheet area (the Nambu-Goto action) 



S[C]^{^T) dt dy^l + n'^T^y^fl^-^j (12) 

generates two types of solutions^ Q 01 ■ One corresponds to a single world-sheet with a nontrivial profile X3{y), 

X3. = ±TTTq r (13) 



where q is a constant of integration determined by the boundary condition 

2q r dy 



^Tjy^ v/(2/4-l)(y'4-y4) 



(14) 



with y'^ — 1 + q^. The corresponding value of ■\/AS'[C] is denoted by Ii. The other solution consists of two parallel 
world sheets with X3 = ±^ extending to the black hole horizon and the corresponding value of •\/A5'[C] is denoted by 
I2 ■ The latter solution corresponds to two non-interacting static quarks in the medium and is equal to the denominator 
of dl 

The free energy we are interested in reads 

F(r, T) = rmin(/, 0) (15) 

where 



/^/i-/2-x/A /°°d2/( ) +l-yc. (16) 



'I 

Inverting ea. (|14p to express r in terms of q and substituting the result to (jl6p . it was found that the function / 
consists of two branches, The upper branch is always positive and is therefore unstable. The lower branch starts 
from being negative for r < tq and becomes positive for r > tq. Both branches joins at r = Tc beyond which the 
nontrivial solution ceases to exist. Numerically, we have ro — "'J^"'" and Vc — Introducing a dimensionless radial 
coordinate, 

p — nTr, (17) 

we find that 

F{r,T) = -%{p)9{po~p), (18) 
r 



where a = ^^^fr^VX ^ 0.2285\/A. The screening factor (/)(p) = — p//(7ra) and is shown in Fig. la. We have 0(0) = 1 

po = 0.7541. (19) 



and 0(po) — with 



The small p expansion of 4>{p) is given by 

p4 (1) op8 (1) 

M-^-^^+^^' + Oi,'; (20, 

On writing the wave function ip{r) ~ ui{p)Yi„i{f) , the radial Schrodinger equation for a zero energy bound state reads 

d^ui 2 dui 



dp^ p dp 



p2 



uz = (21) 
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with V = MKff./(7r2T2). We have 



for the F ansatz and 



pap 



(22) 



V = -- 



Pop 



M-P[^)\0ipo-P) 



for the U ansatz, where 



apoM 



(23) 



(24) 



Note that the potential of the U-ansatz jumps to zero from below at p = po, since the derivative of (/)(p) is nonzero 
there. For both ansatz, and the case with an infrared cutoff discussed below, the solution to ([2T|) is given by [13] 



ui = const. p' 



-i-i 



at p > Po and by 



ui = const, p 



(25) 



(26) 



near the origin. The threshold value of 77 at the dissociation temperature, rjd, is determined by the matching condition 
at p = Po, 



It follows from ((24|) that the dissociation temperature is given by 

It is interesting to observe that the extrapolation of the first two terms of (|20p vanishes at 

47r3 



P^ Po 



0.7178, 



(27) 



(28) 



(29) 



which is very close to the exact zero point (|T9|) . and the third term of ([20]) remains small there. This suggests that 
the screening factor 0(p) can be well approximated by a linear function 



0(p) ^ 1 - ^ 
Po 



(30) 



with Po = ^{po + Po) — 0.7359, as is evident from the exact profile of 0(p) shown in Fig. la. The effective potential 
Ves. is then approximated by a truncated Coulomb potential. We have 



V 



— ( 

pop V 



1 



PO 



0{po - P) 



for the F-ansatz and 



V 



Pop 



-^(Po - P) 



(31) 



(32) 



for the U-ansatz, where the over bar of po has been suppressed. 

The radial wave function of the F-ansatz under the truncated Coulomb approximation can be expressed in terms 
of the confluent hypergeometric function of the 1st kind for p < po, i.e. 



M/=p'iFi(Z + l- J;2/-|-2;2r/-^). 

2 Po 



(33) 
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FIG. 1: (a)The exact screening factor <j>{p) profile extracted from the metric (b) The screening factors r)'s extracted 
from the metric (|47p at different ratios of T/Tc = 1, 2, 3,00 from top to bottom. 



ansatz 


J/^{ls) 


J/^(2s) 


J/^(2p) 


T(ls) 


T(2s) 


T(2p) 


F 


67-124 


15-28 


13-25 


197-364 


44-81 


40-73 


U 


143-265 


27-50 


31-58 


421-779 


80-148 


92-171 



TABLE I: Td in MeV's under the truncated Coulomb approximation. The lower value of each entry corresponds to A = 5.5 
and the upper one to A = 67r. 



The matching condition (|27|) yields the secular equation for 77, 



T] \ iFi{l + 2-^;2l + 3;2r]) 



2i+i-.,+,(^i-^j -;,^._|:„^,:,;; = o. ,34, 

As 77 is reduced from above, we expect the bound state to melt successively. Therefore the first positive root corresponds 
to the minimum binding strength for a bound state of angular momentum / and the 2nd one to the threshold of the 
first radial excitation of the same partial wave. Knowing the values of these Ty's, the disassociation temperature can 
be computed from the formula For example, the threshold ry of the Is state, r/is ~ 1.76, which implies that 

Td ~ 0.0174VAAf. (35) 
In case of the [/-ansatz under the same approximation, we find that 

ui^^J2i+i(2ruP^) (36) 



for p < pq with J,^{x) the Bessel function. The secular equation for 77 reads 

2/ + l-^^^^ = 0. (37) 

We have 771^ — 1.21 and 

Td ~ 0.0368VAAf. (38) 

Numerical results for the dissociation temperature of quarkonium are tabulated in table I, where we have used 
the mass values M — 1.65GeV, 4.85 GeV for c and h quarks [ll[ The errors caused by the truncated Coulomb 
approximation are within 4 percent, as is shown by the numerical solution to the secular equation of the exact 
potential. 

To assess the validity of the potential model employed, we consider a classical two body problem with the truncated 
Coulomb potential (|15p . The Lagrangian of the system reads 

L = ]^M(r\ +>') - Feff.dri - r^]) (39) 

Introducing the center of mass coordinates and the relative coordinates via fi=R + f/2 and r2~R — ^12, and 
assuming a static pair, i? = 0, we have 

L=^Mt^ -V^s.ir) (40) 
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which becomes a particle of the reduced mass in an central potential 14ff.(r). For a circular orbit, the force balance 
at the border of the potential well, r = d = pq/(ttT) is given by 

M ^ ^ (41) 

2d d2 ^ ' 

which yields 

= ^ = 4 (42) 

4 2772 ^ ' 

with —T-f —^2 ^^'^ At melting temperature of U-ansatz, we find that 0.119 < < 0.408 for A in the range 

of ([21). Therefore the approximation may be problematic at the upper limit of the range On the other hand, 
the lower values of a was advocated in [9| based on a comparison between the potentials from lattice simulation and 
that from the AdS/CFT and may serve our purpose better. The NR approximation works better for the dissociation 
temperature of excitations because of the higher 77 values. 



III. QUARKONIUM IN A MODEL WITH AN INFRARED CUTOFF 

Because of the conformal invariance at quantum level, there is no color confinement in = 4 SUSY YM even at 
zero temperature. In order to simulate the confined phase of QCD at low temperature, an infrared cutoff has to be 
introduced that suppress the contribution of the AdS horizon. Two scenarios explored in the literature are the hard- 
wall model and the soft-wall models. The gravity dual of de-confinement transition is modeled as the Hawking-Page 
transition from a metric without a black hole to that with a black hole with T > Tc- The gravity dual of the free 
energy with a hard wall is the Einstein-Hilbert action with a cosmological constant given by 



T 



F = J dzy^iR - 12), (43) 

subject to an appropriate UV regularization, where R is the curvature scalar. In the hadronic phase, the metric 
underlying g and R is that of the standard AdS^ 

ds^ ^ l-(dt^ + dx^ + dz^), (44) 

truncated beyond zq with zq determined by the p- meson mass [l2| . In the plasma phase, the underlying metric is 
given by 

ds^ ^ ^(fdt^ + dx^ + f-^dz^), (45) 
z^ 

with / = 1 — tt'^T^z^, which is identical to ([H]) upon the transformation y = \/{'kTz). The domain of the z-integration 
is < z < 1/{ttT) which corresponds to 1 < y < 00. Notice that zq > 1/{tiT) above the transition temperature 
which was found to be Tc ~ 0.1574mp Therefore heavy quark potential and the meson dissociation temperatures 
calculated above are identical to what calculated above for the vanila A^ = 4 SUSY YM. 

In case of the simplest soft- wall model (jl3l), ^ dilaton is introduced with eq. (|43| modified to 

F = -r-i— / d^x f dre-^^{R-12), (46) 

with c determined by the p-mass and the transition temperature is predicted as Tc ~ 0.2459mp [l^. The string frame 
metric underlying g and R remains given by (|44p in the hadronic phase and by (j45p in the plasma phase. Although, 
the infrared cutoff is partially carried over to the plasma phase, the heavy-quark potential and the dissociation 
temperatures thus obtained remains intact since the minimum area dual to a Wilson loop has to be defined with 
respect to the string frame metric plj . 

A variant of the soft- wall scenario proposed in ref . [isj . however, admits a string frame metric that is different from 
(j45|) by a conformal factor, i. e. 

bz'' 

ds^ = —ifdf + dx^ ^ r^dz^), (47) 
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and 



TTa 



dye^ 



(48) 



The value of 6 = 0.184GeV was obtained by fitting the lattice simulated transition temperature Tc = 186MeV. 
Following the steps from Q to (|16p . we find that 

F{r,T) = -^X{P.T) 
r 

where the screening factor x(p, T) is defined parametrically by 



(49) 



with 



dye 



2q 

■'y- W(y4„ l)(y4_ l_^2g-^) 



2/3 



g^ = (y^-l)e^I 



(50) 



(51) 



and (3 = = 0.539j#-. We have x(p, cx)) = 0(p)- The small p behavior of x(p, T) reads 



x(p,T) = l-^^ap + 0(p2). (52) 

where 

a = e'' - 2v/^ / dxe^ . (53) 
Jo 

The numerical results of x{PiT) for several values of the ratio T/T). are shown in Fig. lb. As is seen, the screening 
becomes weaker in the neighborhood of Tc, similar to lattice QCD [1^. The function x(p, T) vanishes at po — 0.8485 
ai T — Tc while the extrapolation of the first two terms of (|52p vanishes at p'^ ~ 1.764. The truncated Coulomb 
approximation deteriorates in the vicinity of Tc- 

To determine the dissociation temperature in this case, we solve the Schro dinger equation (|2ip numerically with 
the numerically calculated heavy quark potential for both ansatz. We have 

2 

V^-^X{P.T) (54) 
Pop 



for the F ansatz and 

pop 



for the U ansatz. The solution for p < po can be obtained by standard Runge-Kutta method and the threshold value 
of ?7, r]d follows from the matching condition (|27p . Notice that the the infrared cutoff renders both po a-nd r/^ nontrivial 
functions of temperature and eq. lj28p becomes an implicit equation of Td. Nor does Td scales simply with \/A and 
M according to ((35)) and ((38l The modified dissociation temperatures are tabulated in the table II, which show an 
significant increment in the dissociation temperature in the vicinity of T^. The comparison between the ratios T^/Tc 
we calculated here with that obtained from the lattice QCD is shown in table III. 

We observe that the longer screening length near Tc requires that 6 > in the string frame metric (|47p . A negative 
b makes a of (j53p greater than one and Pq smaller than the corresponding value of the conformal ^ . The numerical 
result of po for 6 < is also found to be smaller and the temperature dependence of the screening length is opposite 
to that of the lattice QCD. 
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ansatz 




T 


F 


NA 


235-385 


U 


219-322 


459-779 



TABLE II: Td in MeV's for the Is state with the deformed metric. "NA" means that there is no bound states above Tc and 
the entry for the T with U ansatz and ct = Qn is taken from the table I, since no significant increment is observed. 



ansatz 


J/V' (holographic) 


J/'!/'(lattice) 


T (holographic) 


T (lattice) 


F 


NA 


1.1 


1.3-2.1 


2.3 


U 


1.2-1.7 


2.0 


2.5-4.2 


4.5 



TABLE III: The ratio Td/Tc from the holographic potential (|48j and that from the lattice QCD [T3|. 

IV. CONCLUDING REMARKS 

In summary, we have calculated the dissociation temperatures of heavy quarkonium using the AdS/CFT implied 
potential with both the vanila AdS-Schwarzschild metric and the one with the infrared cutoff. While estimations 
of Td have been made in the literature based on various holographic models [l6l | [l7j . a determination of Td from 
the Schrodinger equation within the same framework is still lacking. Our work is to fill this gap. On comparing 
our results with that from the lattice simulation [IQ] . we found that our ratios T^/Tc extracted from the modified 
AdS-Schwarzschild metric (|47p are lower than the lattice ones within a factor of two. The difference between the 
F-ansatz and the U-ansatz is similar. One has to bear in mind that the matter field proportion of TV = 4 SUSY YM is 
larger than that of QCD with light quarks and the screening ought to be stronger. The stronger screening will make 
the heavy quarkonium more vulnerable and thereby lower the dissociation temperature. This is consistent with the 
observation that the potential becomes wider in the metric with the IR cutoff introduced in [l5| since some of degrees 
of freedom becomes massive. 

The authors of [l^ calculated the spectral function of the fluctuation of a D7 brane and deduced from which the 
meson melting temperature 

T. = ^. (56) 

for all bound state levels. On substituting the value of the heavy quark masses, it gives rise to 825MeV < Td < 1.53GeV 
for J/tp and 2.42GeV < Td < 4.85GeV for T within the range ^ of the 't Hooft coupling. The spectral analysis of 
[18( would be superior if the underlying dynamics of QGP were iV = 4 SUSY YM. But the higher T^'s may point to 
its difference from QCD. Also the A dependence of Td in eq. (|56p is entirely at variance with ours. Our relationship 
Td (X VA for the vanila AdS-Schwarzschild metric follows from the property that the binding strength is proportional 
to -s/A but the binding range is independent of A. A field theoretic speculation on this property of the heavy quark 
potential can be found in [l9|. Since we are comparing two different theories, some features may be shared by both 
and some properties may not. If the screening properties of the = 4 SUSY YM can be carried over to QCD, our 
potential model calculation inspired by the holographic principle should provide a semi-quantitative description on 
the quarkonium dissociation mechanism in the realistic quark-gluon plasma. 
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